Genome sequencing as a generic diagnostic strategy for rare disease.

BACKGROUND: To diagnose the full spectrum of hereditary and congenital diseases, genetic laboratories use many different workflows, ranging from karyotyping to exome sequencing. A single generic high-throughput workflow would greatly increase efficiency. We assessed whether genome sequencing (GS) can replace these existing workflows aimed at germline genetic diagnosis for rare disease. METHODS: We performed short-read GS (NovaSeq™6000; 150 bp paired-end reads, 37 × mean coverage) on 1000 cases with 1271 known clinically relevant variants, identified across different workflows, representative of our tertiary diagnostic centers. Variants were categorized into small variants (single nucleotide variants and indels < 50 bp), large variants (copy number variants and short tandem repeats) and other variants (structural variants and aneuploidies). Variant calling format files were queried per variant, from which workflow-specific true positive rates (TPRs) for detection were determined. A TPR of ≥ 98% was considered the threshold for transition to GS. A GS-first scenario was generated for our laboratory, using diagnostic efficacy and predicted false negative as primary outcome measures. As input, we modeled the diagnostic path for all 24,570 individuals referred in 2022, combining the clinical referral, the transition of the underlying workflow(s) to GS, and the variant type(s) to be detected. RESULTS: Overall, 95% (1206/1271) of variants were detected. Detection rates differed per variant category: small variants in 96% (826/860), large variants in 93% (341/366), and other variants in 87% (39/45). TPRs varied between workflows (79-100%), with 7/10 being replaceable by GS. Models for our laboratory indicate that a GS-first strategy would be feasible for 84.9% of clinical referrals (750/883), translating to 71% of all individuals (17,444/24,570) receiving GS as their primary test. An estimated false negative rate of 0.3% could be expected. CONCLUSIONS: GS can capture clinically relevant germline variants in a 'GS-first strategy' for the majority of clinical indications in a genetics diagnostic lab.

Refining the 9q34.3 microduplication syndrome reveals mild neurodevelopmental features associated with a distinct global DNA methylation profile.

Precise regulation of gene expression is important for correct neurodevelopment. 9q34.3 deletions affecting the EHMT1 gene result in a syndromic neurodevelopmental disorder named Kleefstra syndrome. In contrast, duplications of the 9q34.3 locus encompassing EHMT1 have been suggested to cause developmental disorders, but only limited information has been available. We have identified 15 individuals from 10 unrelated families, with 9q34.3 duplications <1.5 Mb in size, encompassing EHMT1 entirely. Clinical features included mild developmental delay, mild intellectual disability or learning problems, autism spectrum disorder, and behavior problems. The individuals did not consistently display dysmorphic features, congenital anomalies, or growth abnormalities. DNA methylation analysis revealed a weak DNAm profile for the cases with 9q34.3 duplication encompassing EHMT1, which could segregate the majority of the affected cases from controls. This study shows that individuals with 9q34.3 duplications including EHMT1 gene present with mild non-syndromic neurodevelopmental disorders and DNA methylation changes different from Kleefstra syndrome.

Variant-specific pathophysiological mechanisms of AFF3 differently influence transcriptome profiles.

Background: We previously described the KINSSHIP syndrome, an autosomal dominant disorder associated with intellectual disability (ID), mesomelic dysplasia and horseshoe kidney,caused by de novo variants in the degron of AFF3. Mouse knock-ins and overexpression in zebrafish provided evidence for a dominant-negative (DN) mode-of-action, wherein an increased level of AFF3 resulted in pathological effects. Methods: Evolutionary constraints suggest that other mode-of-inheritance could be at play. We challenged this hypothesis by screening ID cohorts for individuals with predicted-to-be deleterious variants in AFF3. We used both animal and cellular models to assess the deleteriousness of the identified variants. Results: We identified an individual with a KINSSHIP-like phenotype carrying a de novo partial duplication of AFF3 further strengthening the hypothesis that an increased level of AFF3 is pathological. We also detected seventeen individuals displaying a milder syndrome with either heterozygous LoF or biallelic missense variants in AFF3. Consistent with semi-dominance, we discovered three patients with homozygous LoF and one compound heterozygote for a LoF and a missense variant, who presented more severe phenotypes than their heterozygous parents. Matching zebrafish knockdowns exhibit neurological defects that could be rescued by expressing human AFF3 mRNA, confirming their association with the ablation of aff3. Conversely, some of the human AFF3 mRNAs carrying missense variants identified in affected individuals did not complement. Overexpression of mutated AFF3 mRNAs in zebrafish embryos produced a significant increase of abnormal larvae compared to wild-type overexpression further demonstrating deleteriousness. To further assess the effect of AFF3 variation, we profiled the transcriptome of fibroblasts from affected individuals and engineered isogenic cells harboring +/+, DN/DN, LoF/+, LoF/LoF or DN/LoF AFF3 genotypes. The expression of more than a third of the AFF3 bound loci is modified in either the DN/DN or the LoF/LoF lines. While the same pathways are affected, only about one-third of the differentially expressed genes are common to these homozygote datasets, indicating that AFF3 LoF and DN variants largely modulate transcriptomes differently, e.g. the DNA repair pathway displayed opposite modulation. Conclusions: Our results and the high pleiotropy shown by variation at this locus suggest that minute changes in AFF3 function are deleterious.

Reanalysis of whole-exome sequencing (WES) data of children with neurodevelopmental disorders in a standard patient care context.

This study aims to inform future genetic reanalysis management by evaluating the yield of whole-exome sequencing (WES) reanalysis in standard patient care in the Netherlands. Single-center data of 159 patients with a neurodevelopmental disorder (NDD), in which WES analysis and reanalysis were performed between January 1, 2014, and December 31, 2021, was retrospectively collected. Patients were included if they were under the age of 18 years at initial analysis and if this initial analysis did not result in a diagnosis. Demographic, phenotypic, and genotypic characteristics of patients were collected and analyzed. The primary outcomes of our study were (i) diagnostic yield at reanalysis, (ii) reasons for detecting a new possibly causal variant at reanalysis, (iii) unsolicited findings, and (iv) factors associated with positive result of reanalysis. In addition, we conducted a questionnaire study amongst the 7 genetic department in the Netherlands creating an overview of used techniques, yield, and organization of WES reanalysis. The single-center data show that in most cases, WES reanalysis was initiated by the clinical geneticist (65%) or treating physician (30%). The mean time between initial WES analysis and reanalysis was 3.7 years. A new (likely) pathogenic variant or VUS with a clear link to the phenotype was found in 20 initially negative cases, resulting in a diagnostic yield of 12.6%. In 75% of these patients, the diagnosis had clinical consequences, as for example, a screening plan for associated signs and symptoms could be devised. Most (32%) of the (likely) causal variants identified at WES reanalysis were discovered due to a newly described gene-disease association. In addition to the 12.6% diagnostic yield based on new diagnoses, reclassification of a variant of uncertain significance found at initial analysis led to a definite diagnosis in three patients. Diagnostic yield was higher in patients with dysmorphic features compared to patients without clear dysmorphic features (yield 27% vs. 6%; p = 0.001). CONCLUSIONS: Our results show that WES reanalysis in patients with NDD in standard patient care leads to a substantial increase in genetic diagnoses. In the majority of newly diagnosed patients, the diagnosis had clinical consequences. Knowledge about the clinical impact of WES reanalysis, clinical characteristics associated with higher yield, and the yield per year after a negative WES in larger clinical cohorts is warranted to inform guidelines for genetic reanalysis. These guidelines will be of great value for pediatricians, pediatric rehabilitation specialists, and pediatric neurologists in daily care of patients with NDD. WHAT IS KNOWN: • Whole exome sequencing can cost-effectively identify a genetic cause of intellectual disability in about 30-40% of patients. • WES reanalysis in a research setting can lead to a definitive diagnosis in 10-20% of previously exome negative cases. WHAT IS NEW: • WES reanalysis in standard patient care resulted in a diagnostic yield of 13% in previously exome negative children with NDD. • The presence of dysmorphic features is associated with an increased diagnostic yield of WES reanalysis.

Mobile element insertions in rare diseases: a comparative benchmark and reanalysis of 60,000 exome samples.

Mobile element insertions (MEIs) are a known cause of genetic disease but have been underexplored due to technical limitations of genetic testing methods. Various bioinformatic tools have been developed to identify MEIs in Next Generation Sequencing data. However, most tools have been developed specifically for genome sequencing (GS) data rather than exome sequencing (ES) data, which remains more widely used for routine diagnostic testing. In this study, we benchmarked six MEI detection tools (ERVcaller, MELT, Mobster, SCRAMble, TEMP2 and xTea) on ES data and on GS data from publicly available genomic samples (HG002, NA12878). For all the tools we evaluated sensitivity and precision of different filtering strategies. Results show that there were substantial differences in tool performance between ES and GS data. MELT performed best with ES data and its combination with SCRAMble increased substantially the detection rate of MEIs. By applying both tools to 10,890 ES samples from Solve-RD and 52,624 samples from Radboudumc we were able to diagnose 10 patients who had remained undiagnosed by conventional ES analysis until now. Our study shows that MELT and SCRAMble can be used reliably to identify clinically relevant MEIs in ES data. This may lead to an additional diagnosis for 1 in 3000 to 4000 patients in routine clinical ES.

DNA methylation episignature, extension of the clinical features, and comparative epigenomic profiling of Hao-Fountain syndrome caused by variants in USP7.

PURPOSE: Hao-Fountain syndrome (HAFOUS) is a neurodevelopmental disorder caused by pathogenic variants in USP7. HAFOUS is characterized by developmental delay, intellectual disability, speech delay, behavioral abnormalities, autism spectrum disorder, seizures, hypogonadism, and mild dysmorphic features. We investigated the phenotype of 18 participants with HAFOUS and performed DNA methylation (DNAm) analysis, aiming to generate a diagnostic biomarker. Furthermore, we performed comparative analysis with known episignatures to gain more insight into the molecular pathophysiology of HAFOUS. METHODS: We assessed genomic DNAm profiles of 18 individuals with pathogenic variants and variants of uncertain significance (VUS) in USP7 to map and validate a specific episignature. The comparison between the USP7 cohort and 56 rare genetic disorders with earlier reported DNAm episignatures was performed with statistical and functional correlation. RESULTS: We mapped a sensitive and specific DNAm episignature for pathogenic variants in USP7 and utilized this to reclassify the VUS. Comparative epigenomic analysis showed evidence of HAFOUS similarity to a number of other rare genetic episignature disorders. CONCLUSION: We discovered a sensitive and specific DNAm episignature as a robust diagnostic biomarker for HAFOUS that enables VUS reclassification in USP7. We also expand the phenotypic spectrum of 9 new and 5 previously reported individuals with HAFOUS.

PhenoScore quantifies phenotypic variation for rare genetic diseases by combining facial analysis with other clinical features using a machine-learning framework.

Several molecular and phenotypic algorithms exist that establish genotype-phenotype correlations, including facial recognition tools. However, no unified framework that investigates both facial data and other phenotypic data directly from individuals exists. We developed PhenoScore: an open-source, artificial intelligence-based phenomics framework, combining facial recognition technology with Human Phenotype Ontology data analysis to quantify phenotypic similarity. Here we show PhenoScore's ability to recognize distinct phenotypic entities by establishing recognizable phenotypes for 37 of 40 investigated syndromes against clinical features observed in individuals with other neurodevelopmental disorders and show it is an improvement on existing approaches. PhenoScore provides predictions for individuals with variants of unknown significance and enables sophisticated genotype-phenotype studies by testing hypotheses on possible phenotypic (sub)groups. PhenoScore confirmed previously known phenotypic subgroups caused by variants in the same gene for SATB1, SETBP1 and DEAF1 and provides objective clinical evidence for two distinct ADNP-related phenotypes, already established functionally.

Comprehensive de novo mutation discovery with HiFi long-read sequencing.

BACKGROUND: Long-read sequencing (LRS) techniques have been very successful in identifying structural variants (SVs). However, the high error rate of LRS made the detection of small variants (substitutions and short indels < 20 bp) more challenging. The introduction of PacBio HiFi sequencing makes LRS also suited for detecting small variation. Here we evaluate the ability of HiFi reads to detect de novo mutations (DNMs) of all types, which are technically challenging variant types and a major cause of sporadic, severe, early-onset disease. METHODS: We sequenced the genomes of eight parent-child trios using high coverage PacBio HiFi LRS (~ 30-fold coverage) and Illumina short-read sequencing (SRS) (~ 50-fold coverage). De novo substitutions, small indels, short tandem repeats (STRs) and SVs were called in both datasets and compared to each other to assess the accuracy of HiFi LRS. In addition, we determined the parent-of-origin of the small DNMs using phasing. RESULTS: We identified a total of 672 and 859 de novo substitutions/indels, 28 and 126 de novo STRs, and 24 and 1 de novo SVs in LRS and SRS respectively. For the small variants, there was a 92 and 85% concordance between the platforms. For the STRs and SVs, the concordance was 3.6 and 0.8%, and 4 and 100% respectively. We successfully validated 27/54 LRS-unique small variants, of which 11 (41%) were confirmed as true de novo events. For the SRS-unique small variants, we validated 42/133 DNMs and 8 (19%) were confirmed as true de novo event. Validation of 18 LRS-unique de novo STR calls confirmed none of the repeat expansions as true DNM. Confirmation of the 23 LRS-unique SVs was possible for 19 candidate SVs of which 10 (52.6%) were true de novo events. Furthermore, we were able to assign 96% of DNMs to their parental allele with LRS data, as opposed to just 20% with SRS data. CONCLUSIONS: HiFi LRS can now produce the most comprehensive variant dataset obtainable by a single technology in a single laboratory, allowing accurate calling of substitutions, indels, STRs and SVs. The accuracy even allows sensitive calling of DNMs on all variant levels, and also allows for phasing, which helps to distinguish true positive from false positive DNMs.

A complex structural variant near SOX3 causes X-linked split-hand/foot malformation.

Split-hand/foot malformation (SHFM) is a congenital limb defect most typically presenting with median clefts in hands and/or feet, that can occur in a syndromic context as well as in isolated form. SHFM is caused by failure to maintain normal apical ectodermal ridge function during limb development. Although several genes and contiguous gene syndromes are implicated in the monogenic etiology of isolated SHFM, the disorder remains genetically unexplained for many families and associated genetic loci. We describe a family with isolated X-linked SHFM, for which the causative variant could be detected after a diagnostic journey of 20 years. We combined well-established approaches including microarray-based copy number variant analysis and fluorescence in situ hybridization coupled with optical genome mapping and whole genome sequencing. This strategy identified a complex structural variant (SV) comprising a 165-kb gain of 15q26.3 material ([GRCh37/hg19] chr15:99795320-99960362dup) inserted in inverted position at the site of a 38-kb deletion on Xq27.1 ([GRCh37/hg19] chrX:139481061-139518989del). In silico analysis suggested that the SV disrupts the regulatory framework on the X chromosome and may lead to SOX3 misexpression. We hypothesize that SOX3 dysregulation in the developing limb disturbed the fine balance between morphogens required for maintaining AER function, resulting in SHFM in this family.

Rapid exome sequencing as a first-tier test in neonates with suspected genetic disorder: results of a prospective multicenter clinical utility study in the Netherlands.

The introduction of rapid exome sequencing (rES) for critically ill neonates admitted to the neonatal intensive care unit has made it possible to impact clinical decision-making. Unbiased prospective studies to quantify the impact of rES over routine genetic testing are, however, scarce. We performed a clinical utility study to compare rES to conventional genetic diagnostic workup for critically ill neonates with suspected genetic disorders. In a multicenter prospective parallel cohort study involving five Dutch NICUs, we performed rES in parallel to routine genetic testing for 60 neonates with a suspected genetic disorder and monitored diagnostic yield and the time to diagnosis. To assess the economic impact of rES, healthcare resource use was collected for all neonates. rES detected more conclusive genetic diagnoses than routine genetic testing (20% vs. 10%, respectively), in a significantly shorter time to diagnosis (15 days (95% CI 10-20) vs. 59 days (95% CI 23-98, p < 0.001)). Moreover, rES reduced genetic diagnostic costs by 1.5% (€85 per neonate). CONCLUSION:  Our findings demonstrate the clinical utility of rES for critically ill neonates based on increased diagnostic yield, shorter time to diagnosis, and net healthcare savings. Our observations warrant the widespread implementation of rES as first-tier genetic test in critically ill neonates with disorders of suspected genetic origin. WHAT IS KNOWN: • Rapid exome sequencing (rES) enables diagnosing rare genetic disorders in a fast and reliable manner, but retrospective studies with neonates admitted to the neonatal intensive care unit (NICU) indicated that genetic disorders are likely underdiagnosed as rES is not routinely used. • Scenario modeling for implementation of rES for neonates with presumed genetic disorders indicated an expected increase in costs associated with genetic testing. WHAT IS NEW: • This unique prospective national clinical utility study of rES in a NICU setting shows that rES obtained more and faster diagnoses than conventional genetic tests. • Implementation of rES as replacement for all other genetic tests does not increase healthcare costs but in fact leads to a reduction in healthcare costs.

Whole genome sequencing for USH2A-associated disease reveals several pathogenic deep-intronic variants that are amenable to splice correction.

A significant number of individuals with a rare disorder such as Usher syndrome (USH) and (non-)syndromic autosomal recessive retinitis pigmentosa (arRP) remain genetically unexplained. Therefore, we assessed subjects suspected of USH2A-associated disease and no or mono-allelic USH2A variants using whole genome sequencing (WGS) followed by an improved pipeline for variant interpretation to provide a conclusive diagnosis. One hundred subjects were screened using WGS to identify causative variants in USH2A or other USH/arRP-associated genes. In addition to the existing variant interpretation pipeline, a particular focus was put on assessing splice-affecting properties of variants, both in silico and in vitro. Also structural variants were extensively addressed. For variants resulting in pseudoexon inclusion, we designed and evaluated antisense oligonucleotides (AONs) using minigene splice assays and patient-derived photoreceptor precursor cells. Biallelic variants were identified in 49 of 100 subjects, including novel splice-affecting variants and structural variants, in USH2A or arRP/USH-associated genes. Thirteen variants were shown to affect USH2A pre-mRNA splicing, including four deep-intronic USH2A variants resulting in pseudoexon inclusion, which could be corrected upon AON treatment. We have shown that WGS, combined with a thorough variant interpretation pipeline focused on assessing pre-mRNA splicing defects and structural variants, is a powerful method to provide subjects with a rare genetic condition, a (likely) conclusive genetic diagnosis. This is essential for the development of future personalized treatments and for patients to be eligible for such treatments.

The Genetics of Intellectual Disability.

Intellectual disability (ID) has a prevalence of ~2-3% in the general population, having a large societal impact. The underlying cause of ID is largely of genetic origin; however, identifying this genetic cause has in the past often led to long diagnostic Odysseys. Over the past decades, improvements in genetic diagnostic technologies and strategies have led to these causes being more and more detectable: from cytogenetic analysis in 1959, we moved in the first decade of the 21st century from genomic microarrays with a diagnostic yield of ~20% to next-generation sequencing platforms with a yield of up to 60%. In this review, we discuss these various developments, as well as their associated challenges and implications for the field of ID, which highlight the revolutionizing shift in clinical practice from a phenotype-first into genotype-first approach.

A Solve-RD ClinVar-based reanalysis of 1522 index cases from ERN-ITHACA reveals common pitfalls and misinterpretations in exome sequencing.

PURPOSE: Within the Solve-RD project (https://solve-rd.eu/), the European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies aimed to investigate whether a reanalysis of exomes from unsolved cases based on ClinVar annotations could establish additional diagnoses. We present the results of the "ClinVar low-hanging fruit" reanalysis, reasons for the failure of previous analyses, and lessons learned. METHODS: Data from the first 3576 exomes (1522 probands and 2054 relatives) collected from European Reference Network for Intellectual disability, TeleHealth, Autism and Congenital Anomalies was reanalyzed by the Solve-RD consortium by evaluating for the presence of single-nucleotide variant, and small insertions and deletions already reported as (likely) pathogenic in ClinVar. Variants were filtered according to frequency, genotype, and mode of inheritance and reinterpreted. RESULTS: We identified causal variants in 59 cases (3.9%), 50 of them also raised by other approaches and 9 leading to new diagnoses, highlighting interpretation challenges: variants in genes not known to be involved in human disease at the time of the first analysis, misleading genotypes, or variants undetected by local pipelines (variants in off-target regions, low quality filters, low allelic balance, or high frequency). CONCLUSION: The "ClinVar low-hanging fruit" analysis represents an effective, fast, and easy approach to recover causal variants from exome sequencing data, herewith contributing to the reduction of the diagnostic deadlock.

The performance of genome sequencing as a first-tier test for neurodevelopmental disorders.

Genome sequencing (GS) can identify novel diagnoses for patients who remain undiagnosed after routine diagnostic procedures. We tested whether GS is a better first-tier genetic diagnostic test than current standard of care (SOC) by assessing the technical and clinical validity of GS for patients with neurodevelopmental disorders (NDD). We performed both GS and exome sequencing in 150 consecutive NDD patient-parent trios. The primary outcome was diagnostic yield, calculated from disease-causing variants affecting exonic sequence of known NDD genes. GS (30%, n = 45) and SOC (28.7%, n = 43) had similar diagnostic yield. All 43 conclusive diagnoses obtained with SOC testing were also identified by GS. SOC, however, required integration of multiple test results to obtain these diagnoses. GS yielded two more conclusive diagnoses, and four more possible diagnoses than ES-based SOC (35 vs. 31). Interestingly, these six variants detected only by GS were copy number variants (CNVs). Our data demonstrate the technical and clinical validity of GS to serve as routine first-tier genetic test for patients with NDD. Although the additional diagnostic yield from GS is limited, GS comprehensively identified all variants in a single experiment, suggesting that GS constitutes a more efficient genetic diagnostic workflow.

Optical genome mapping and revisiting short-read genome sequencing data reveal previously overlooked structural variants disrupting retinal disease-associated genes.

PURPOSE: Structural variants (SVs) play an important role in inherited retinal diseases (IRD). Although the identification of SVs significantly improved upon the availability of genome sequencing, it is expected that involvement of SVs in IRDs is higher than anticipated. We revisited short-read genome sequencing data to enhance the identification of gene-disruptive SVs. METHODS: Optical genome mapping was performed to improve SV detection in short-read genome sequencing-negative cases. In addition, reanalysis of short-read genome sequencing data was performed to improve the interpretation of SVs and to re-establish SV prioritization criteria. RESULTS: In a monoallelic USH2A case, optical genome mapping identified a pericentric inversion (173 megabase), with 1 breakpoint disrupting USH2A. Retrospectively, the variant could be observed in genome sequencing data but was previously deemed false positive. Reanalysis of short-read genome sequencing data (427 IRD cases) was performed which yielded 30 pathogenic SVs affecting, among other genes, USH2A (n = 15), PRPF31 (n = 3), and EYS (n = 2). Eight of these (>25%) were overlooked during previous analyses. CONCLUSION: Critical evaluation of our findings allowed us to re-establish and improve our SV prioritization and interpretation guidelines, which will prevent missing pathogenic events in future analyses. Our data suggest that more attention should be paid to SV interpretation and the current contribution of SVs in IRDs is still underestimated.

De novo mutation hotspots in homologous protein domains identify function-altering mutations in neurodevelopmental disorders.

Variant interpretation remains a major challenge in medical genetics. We developed Meta-Domain HotSpot (MDHS) to identify mutational hotspots across homologous protein domains. We applied MDHS to a dataset of 45,221 de novo mutations (DNMs) from 31,058 individuals with neurodevelopmental disorders (NDDs) and identified three significantly enriched missense DNM hotspots in the ion transport protein domain family (PF00520). The 37 unique missense DNMs that drive enrichment affect 25 genes, 19 of which were previously associated with NDDs. 3D protein structure modeling supports the hypothesis of function-altering effects of these mutations. Hotspot genes have a unique expression pattern in tissue, and we used this pattern alongside in silico predictors and population constraint information to identify candidate NDD-associated genes. We also propose a lenient version of our method, which identifies 32 hotspot positions across 16 different protein domains. These positions are enriched for likely pathogenic variation in clinical databases and DNMs in other genetic disorders.

Lessons learned from rapid exome sequencing for 575 critically ill patients across the broad spectrum of rare disease.

Introduction: Rapid exome sequencing (rES) has become the first-choice genetic test for critically ill patients, mostly neonates, young infants, or fetuses in prenatal care, in time-sensitive situations and when it is expected that the genetic test result may guide clinical decision making. The implementation of rES has revolutionized medicine by enabling timely identification of genetic causes for various rare diseases. The utilization of rES has increasingly been recognized as an essential diagnostic tool for the identification of complex and undiagnosed genetic disorders. Methods: We conducted a retrospective evaluation of our experiences with rES performed on 575 critically ill patients from various age groups (prenatal to adulthood), over a four-year period (2016-2019). These patients presented with a wide spectrum of rare diseases, including but not limited to neurological disorders, severe combined immune deficiency, and cancer. Results: During the study period, there was a significant increase in rES referrals, with a rise from a total of two referrals in Q1-2016 to 10 referrals per week in Q4-2019. The median turnaround time (TAT) decreased from 17 to 11 days in the period 2016-2019, with an overall median TAT of 11 days (IQR 8-15 days). The overall diagnostic yield for this cohort was 30.4%, and did not significantly differ between the different age groups (e.g. adults 22.2% vs children 31.0%; p-value 0.35). However, variability in yield was observed between clinical entities: craniofacial anomalies yielded 58.3%, while for three clinical entities (severe combined immune deficiency, aneurysm, and hypogonadotropic hypogonadism) no diagnoses were obtained. Discussion: Importantly, whereas clinical significance is often only attributed to a conclusive diagnosis, we also observed impact on clinical decision-making for individuals in whom no genetic diagnosis was established. Hence, our experience shows that rES has an important role for patients of all ages and across the broad spectrum of rare diseases to impact clinical outcomes.

Episignature Mapping of TRIP12 Provides Functional Insight into Clark-Baraitser Syndrome.

Clark-Baraitser syndrome is a rare autosomal dominant intellectual disability syndrome caused by pathogenic variants in the TRIP12 (Thyroid Hormone Receptor Interactor 12) gene. TRIP12 encodes an E3 ligase in the ubiquitin pathway. The ubiquitin pathway includes activating E1, conjugating E2 and ligating E3 enzymes which regulate the breakdown and sorting of proteins. This enzymatic pathway is crucial for physiological processes. A significant proportion of TRIP12 variants are currently classified as variants of unknown significance (VUS). Episignatures have been shown to represent a powerful diagnostic tool to resolve inconclusive genetic findings for Mendelian disorders and to re-classify VUSs. Here, we show the results of DNA methylation episignature analysis in 32 individuals with pathogenic, likely pathogenic and VUS variants in TRIP12. We identified a specific and sensitive DNA methylation (DNAm) episignature associated with pathogenic TRIP12 variants, establishing its utility as a clinical biomarker for Clark-Baraitser syndrome. In addition, we performed analysis of differentially methylated regions as well as functional correlation of the TRIP12 genome-wide methylation profile with the profiles of 56 additional neurodevelopmental disorders.